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Abstract 

Some optical experiments provide the easiest way to test quantum mechanical predictions. Such 
a situation applies to a laser beam traversing a dielectric structure. The dielectric structure mir- 
roring quantum mechanical potentials. The simplest of these, studied in this paper, involves a 
single dielectric block. We exhibit both by analytic and numerical calculations explicit examples 
of total coherence and total incoherence phenomena. Which case appears depends upon the angle 
of incidence of the laser beam and/or the dimensions of beam vs block width. Unlike our previ- 
ous studies our calculations here are three-dimensional, although, with suitable approximations, 
somewhat simpler planar expressions can be derived. 

I. INTRODUCTION 

Everything in nature is quantum mechanical at its core but, as the history of Physics teaches us, many 
processes can be studied through "classical" equations. More specifically vifith equations that do not 
involve explicitly (or implicitly as occurs when using natural units) Planck's constant. Examples 
are Newton's equation, thermodynamic equations, classical field equations and amongst these, of 
particular interest to us. Maxwell's equations [1]. Maxwell's equations encompass optical phenomena 
which have a surprising aspect. Under certain conditions [2,3], which we recall in the next section, 
the equations of say the electric field is analogous to the non relativistic Schrodinger equation [4]. 
The potentials in Quantum Mechanics (QM) are identifiable with dielectric materials with different 
refractive indices with which the laser interacts. One of the advantages of this fact is that QM effects, 
probabilistic in nature and valid for massive particle wave packets, can be reproduced, and the theory 
tested far easier, with optical analog experiments [5-12]. 

In this paper we wish to simulate one of the simplest of all QM potentials, i.e. a square potential. 
We present both analytic and numerical calculations of a gaussian laser incident upon a single block 
dielectric. We shall in particular concentrate upon the transmitted beam or beams. It will be shown 
that in the so-called wave limit, where coherence between individual amplitude terms dominates, a 
single outgoing laser beam emerges. In the other limit, when incoherence dominates, multiple laser 
beams are transmitted. We shall calculate both the intensity and position of the first few of these 
beams. Our intention is to apply these methods to more complex systems, such as twin barrier 
diffusion and tunnelling [13,14]. The basic ingredients of our approach will be detailed in this paper. 
This work differs from some previous papers of the present authors [15-18] by calculating in three 
dimensions and by the explicit use of a laser as the light source. This provides for a more realistic 
theoretical study, closer to laboratory experiments. For some of the analytical calculations we shall 
make use of the stationary phase method, SPM [19,20]. We shall demonstrate in detail how our 
results depend critically upon the laser angle of incidence with the dielectric block. 

In the next section we derive the gaussian laser amplitude in three-dimensional free space. In 
section III we describe the geometry of the laser with respect to a dielectric block corresponding to a 
QM well. Also in this section we derive the modification to the laser beam within the dielectric. In 
section IV we derive the transmission and reflection amplitudes for the laser beam (or beams) exiting 
the dielectric. The method of calculation follows the so called two step approach [15-17]. This yields 
an infinite number of contributing amplitudes. We study the conditions for which complete coherence 
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occurs yielding a single outgoing beam. In the other extreme infinite outgoing beams are predicted. 
We calculate the exit points and intensities of the leading contributions as a function of the incoming 
angle. In section V we calculate the first few incoherent beam intensities and their exit points from 
the dielectric block. This is performed both analytically and numerically. For this analysis typical 
helium- neon laser parameters are employed [21]. 



II. THE GAUSSIAN LASER AMPLITUDE 



Throughout this paper we consider a fixed frequency w = |fc|/c. A fixed frequency implies a stationary 
system, so we can forego writing the time dependence e~*"*. For the theoretical analysis we start 
with plane waves. These correspond to momentum eigenvalues and determine a fixed direction k. A 
general solution can then be obtained by convoluting plane waves with a suitable integral. 

Consider a plane wave moving in the fc-dircction, the electric field vector can be then expressed 

by 

where s is a unitary polarization vector, assumed linear, which lies in the plane perpendicular to k. 
The condition div£^ = imposed by one of the Maxwell equations implies k ■ § = 0. This means that 
s is k dependent. Now a possible but not unique choice for s is 



~k 



k. 



We claim this is not unique since, as we vary the direction, we may also rotate the above s about 
the plane wave direction with impunity. Furthermore, as given, these do not determine a unique s in 



the limit k^. 



0. In fact the direction of s depends upon how these limits are performed. If first 



we take kx ^ and then fcj, — )■ we find s = (0, 1, 0). If we first take fcj, — J- and then fc^; — >■ we 
find s = (1,0,0). A general direction in the xy-plane is obtained by fixing the ratio of kx/ky and 
taking simultaneously k^^y 0. All of the above tells us that the components of E are not the same 
and thus, after integrations with a chosen convolution function, they will yield different component 
functions of E. However, if use is made, as we shall do, of a peaked gaussian function that limits all 
significant values of wave number to k^^y <^ kz ~ k we can approximate s by a fixed vector in the 
xy-plane. The calculation of E then greatly simplifies to the calculation of the amplitude E. 

In free space the gaussian laser is obtained by convoluting a plane wave with a peaked gaussian, 

G{kx, ky) = exp 

where the arguments are the components of the wave number in the direction orthogonal to the laser 
direction, here chosen as the z direction, and with a Dirac delta function. 
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S{kz - ^k- ~ k-^- ky) . 
The electric field amplitude, E{x,y,z), is then proportional to 



d'''fc G{kx,ky) S{kz — \J k^ — kl + ky ) exp [i {k^ x + ky y + k^ z)] . 

The integral over dfc^ is straightforward because of the Dirac delta function. So we have 

w" f 

E{x,y,z) = Eo^ dkxdky G{kx,ky) exp[i {kxX + kyy + kz z)] , 



(1) 



where Eq = E(0,0,0). The use of a sharp gaussian also allows us to approximate the electric field 
amplitude above as follows 



W" fc f 

E{x, y, z)^ Eo-^ e**^^' / dfc^ dky G{kx,ky) exp 
47r J 



i \ kx X -\- ky y 



(2) 



From this, after performing two generalized gaussian integrations, we obtain the weU known closed 
formula which describes the free propagation of a gaussian laser beam [22], i.e. 



E{x, y, z) sa exp[ ik z] 



w^ + 2 i z/ A: 



exp 



X +y 



wl + 2i z/ k 



(3) 



The intensity distribution of the gaussian beam, /(x, y, z) = \E{x, y,x)\ , is then given by 



I{x,y,z) 



2P 



TT W^(z) 



exp 



X +y 

w^{z) 



(4) 



where 



W(Z) Wo W 1 + 



is the radius of the 1 /e^ contour after the wave has propagated a distance z and P the total power in 
the beam, 

P = y dxdyI{x,y,z)^^Io , 

with /o = /(0,0,0). 

In Fig 1 we plot a side view of the laser intensity for different beam waists, Wq = 1 mm and 
Wo = 2 mm, and for a wavelength value of A = 27r/fc = 633 nm. These figures display the fact that the 
spread in momentum is complementary to the spatial laser width. Indeed with the angular opening 
9 defined as 

w(z) A 



tan^^ = lim 



TT Wo 



the 1/e^ intensity contours asymptotically approach a cone of angular radius 9. For our choice of 
lambda and for a beam waist of the order of mm (A <^ nwo) the angular radius is then given by 

9 X / TTWo . 

This angle is indicative of the spread in wavelength and as given by G{kj., ky) decreases with increasing 
Wq. On the contrary, as seen by the explicit form of the intensity I{x, y, z), the width of a laser, for any 
chosen value of z, grows with Wq. As an aside we recall that the wave number and position variance are 
not "complementary" (in the QM sense) since there is, of course, no uncertainty principle involved for 
classical variables. Finally we observe that after convolution the laser rays do not travel in a unique 
direction. However we shall refer to the laser direction as the direction of the ray of peak intensity. 
This is the z-axis in Eq. (3) . 



III. GEOMETRY OF A PROPOSED EXPERIMENT 



We sketch below the geometry of a laser source at O, with the z-axis extended to a dielectric block 
with normal the z^^,-axls. The angle between the z and z^^ axis is a. A laser impinging along the z-axis 
travels, initially, within the dielectric with refractive index greater than air in a different direction, 
Zd, with angle /3(< a) w.r.t the z*-axis and angle a — (3 w.r.t. the z-axis. The outgoing situation 
when the laser exits the dielectric is, as we shall see, conditioned by the impact angle and the laser 
spatial dimensions. The axis indicated by and z*, and for completeness x*, are obtained from y 
and z, and x, by a rotation of a around the x-axis. We chose the z-axis as that of the center of the 
laser beam so that we may use directly the expressions for the laser of the previous section. Since the 
dielectric block is perpendicular to the z*-axis, it is stratified in the z* direction. 



Air 




If a plane wave traverses the block, the separation of variables implies that the wave numbers in 
the a;* (= x) and directions remain unaltered. Only the wave number changes, in analogy to the 
momentum in one-dimensional QM, thus within the dielectric we have the following plane wave 

Such a plane wave satisfies the Maxwell equations for a stratified media, 

(c*!. +C +C) ^ + "'(^*)^'£^==0 . (5) 

From this equation we get 

Note that for n > 1, is a real number. Thus, in our setup, for whatever incident angle there 
are no evanescent waves and hence no tunneling. We have only diffusion in QM terms. This implies 
that our dielectric block corresponds to a potential well. Now we observe that, after convolution with 

G[kx, ky), 

kz, = kz cos a ~ ky sin a ~ k cos a — ky sin a , 

where the approximation is justified if we neglect terms of second order in k^ and ky and consequently 
spreading effects. Similarly we find, after using Snell's law sin a = n sin/3, 

Qzt ~ nk cos P — ky tan /3 cos a . (6) 

Using these approximations and the fact that k r ~ k^ r^^ the plane wave phase within the dielectric, 

kx^ ^* ~t- ky^ -\- (jz^ z^ ? 

can be rewritten as follows 

kx X -\- ky y -\- kz z -\- {qz, kz^ ) 2* 

sa kx X + ky [y + (sin a — tan /? cos a) ] + fc [z + {n cos /3 — cos a) ] 

= kxX + ky [ cos(a - ^) y + sin(a - I3)z]+nk [ cos(a - (3) z - sin(a - (i)y] ■ 

Vd 2ti 

However we should not simply convolute the above plane wave with G{kx,ky). Wc have first to 
calculate the '"step" refiection and transmission amplitudes. Indeed there is a reflected beam at the 



interface. Nevertlieless, if tlie transmission amplitude at the interface is close to one we can ignore 
reflection and write 



Edix,y,z) E„ exp[inkzd] exp 



X + (cos a/ cos 



(7) 



As expected this is seen to represent a laser with principal direction along the z^-axis and with the 
transverse yd spread increased by cos/3/ cos a. This deforms the laser cross section from circular to 
elliptic. This result can also be directly derived by simple geometry and Snell's law. Our analytic 
derivation confirms the validity of our approximations which will also be used in later sections. 

In the next section we calculate the reflection and transmission coefficients with the method in 
which entry and exit from the dielectric are treated separately as successive and repeated "step" 
calculations [15]. 



IV. TRANSMISSION AND REFLECTION AMPLITUDES 



As we have demonstrated in previous works [16-18], the transmission and reflection amplitudes both 
in Optics and QM can be built up by successive applications of the "step" analysis at dielectric or 
potential discontinuities. In this case at dielectric interfaces. The general z^, solutions, be they oscil- 
latory or evanescent, are constrained by applying continuity to the amplitude and its first derivative 
in z». 

This procedure is repeated ad infinitum for multiply reflected amplitudes. The final "formal" 
transmission and reflected amplitudes are obtained by summing the appropriate series. However if, 
for example, incoherence dominates between successive contributions then the individual terms are 
more relevant than the formal sum. In numerical calculations, which involve the convolution function 
and the plane wave as well, we must work, in general, with the summed amplitudes to reproduce any, 
even partial, coherence phenomena. Our two-step calculations are one-dimcnsional, z*, and perfectly 
recall those of QM. 




The above trilogy of figures are those relevant to a plane wave impinging upon a dielectric block. In 
the first figure (step 1 — > n) i?i is the reflected amplitude and the transmitted amplitude within 
the dielectric for an incident plane wave of unit amplitude. The continuity equations when solved 
yield 

i?i = cxp \ 2ikz, a] , = 2 — exp [ i (kz, — qz ) a] . (8) 

kz, + qz, kz, + qz. 

Note that, in accordance with flux conservation [1], 



The second figure (step n — >■ 1) corresponds to the exiting of the plane wave from the dielectric. Again 
it is convenient to use a unit normalization for the incoming wave. The results for the reflected, R2, 
and transmitted, T2, amplitudes are given by the previous expressions with kz_^ O Qz, and a — > a + 6, 

R2 ^ I"' ~ cxp[2zg,, (a + b)] , T, = 2—^ exp [ i (q,. - fc, J (a + 6) ] . (9) 

Again we observe that, as required, 

\R.\ + — |r,| =1 . 

Finally a reflected wave at the second interface, = a + &, returns to the first interface, = a, and 
again gives rise to a reflected, and transmitted, Ti, wave. The solution of the continuity equations 
are now obtained from i?i and Ti by kz, ~ Qz^, 

-Ri = 7^^— — ^ exp [ - 2 i q^, a ] , fi=2--—^ cxp [i (kz, - qz,) a] . (10) 

These amplitudes satisfy 

, ~ ,2 k, , ~ ,2 

l^il + =1 . 

HZ, 

The total transmission and reflection amplitudes can now be built up by repeated multiplication and 
summation of the reflection and transmission coefficients just given. The figure below gives the first 
few of these multiple reflection and transmission terms. 




The first contribution to the total reflection amplitude is simply Ri . The second contribution is clearly 
T1R2T1. The third is T1R2R1R2T1 and so forth. Hence the full reflection amplitude is 



MT^.- T1R2T1 i{q^^ - k^J sin{qz,b)e'^ 

R = Ri + T1R2 y {R1R2) = Ri ^ = — — — —2 — ; — 2 s ■ / , ^ 

1 - i?ii?2 ^kz,qz, COs{qz,b) - i [k^^ + q^Jsm{qz,b) 

In the same way, it is easy to derive the full transmitted amplitude 

X = T (R R ) T = Z '^ = '^z,qz, e 

' ' ' I-R1R2 '2kz,qz,cos{qz,b)-i{kl^ + qljsin{qz,b) ' 

Now these amplitudes, after some calculation, satisfy the well known condition [1,4] 

\Rf + \Tf = 1 , 



(11) 



(12) 



as must be the case when total coherence and a single reflected and transmitted wave in free space 
occur. However we also find that the sum of the modulus squared of all the contributions to reflection, 

,2 ^ \TMf _ jk,, - q,S 

i~\RiR2\' K + qI 

and to transmission, 

1- |i?ii?,r K, + ql ' 

also sum to one. This result is reasonable since it implies conservation of intensity in the limit of 
total incoherence. Obviously a similar result must hold for any level of coherence. We recall that 
the transmission and reflection intensities are crucially dependent upon the degree of coherence. For 
example, total coherence displays the phenomena of resonance [4]. 



V. TRANSMISSION PHENOMENA 



In this section we calculate both numerically and with the use of the SPM the exit positions of 
individual beams for various incident angles of our gaussian laser. For these calculations, we have 
chosen 

A ^ 633 nm , w^ = 2 mm , 

a = 10 cm , b = 5 cm , 

n = V3 , a = 7r/60 ,7r/30,7r/6,7r/3 . 

The expression for the transmitted amplitude is, with the approximation described in the text for k^, 



ET.{x,y,z) = Ea 



47r 



dkx dky T{kx, ky) G{kx,ky) exp [i {k^ x + ky y + k^ z)] 



(13) 



Now, as shown in Table 1, the dependence of T upon k^ is almost negligible. This follows from the fact 
that all angles are rotations about the x-axis and thus do not involve k^ ■ On the other hand the angles 
do depend upon ky to first order. Thus we can approximate T{kx,ky) with T{0,ky). Furthermore, 
the use of a sharp gaussian allows us to use k^ ~ k — (k^. + ky)/2k. Consequently the k^ integration 
can be performed analytically, as for the incoming laser, and we obtain 



ET,{x,y,z) 



E„ ^ e'''' 



2iz/ k 



exp 



dky r(0, ky) exp 



2k 



2iz/ k 
\ + ikyy 



(14) 



This integral must be performed numerically. In Fig. 2, we plot the results for \Et/Eo\ versus y for 
four values of the incoming a. The x-dependence is not plotted because it is the same as the incoming 
laser beam. We observe that for small a angles a single beam emerges. However, its peak is somewhat 
displaced w.r.t. the incoming beam (an effect of T). As a increases, evidence of a secondary beam 
is seen. For even high a, this beam and others separate and we observe the effect of the incoherent 
limit of the transmission coefficient. Our curves are calculated directly from the full T, however they 
of course agree with those calculated from each appropriate transmission term. We exhibit only the 
first three terms in our plots. 

In Fig. 3 we plot these same results in a different manner suitable to direct experimental verifica- 
tion. We plot the xy-profdc of each intensity for the four chosen a. For these plots we have arbitrarily 
chosen z = 1 m, however, with the very small spreading in our simulation, the exact value of z (below 
several meters) is irrelevant. All centers lie upon the x-axis as does the incoming beam profile. How- 
ever the y values depend upon a. Another way of determining these y values without going through 
a numerical calculation is by means of the SPM. We must now work with the individual transmission 
terms since the SPM only yields a single peak position y. The SPM applied to an expression which 



involves a sum of peaks yields only the mean value, normally not even coincident with any of the 
peaks. The SPM involves taking the derivative of the argument of the phase of a peaked function 
assuming it more relevant to the peak's position compared to the modulus of the amplitude. Now the 
phase of our amplitudes take the form 

^ilkyy+f{ky)] 

where f{ky) is a function of ky (independent of any y dependence). The function / is calculated 
below. The SPM then says that the peak occurs where 



dfiky) 



dk,. 



. 



Consider, the first leading transmission amplitude term T1T2, the phase is 

{qz, ~ k^_^){a + b) + [k^^ - qzja= [q^^ - k^Jb , 

to which, we must add the plane wave phase e'*'"^. Thus in this case, f{ky) — {qz, — k^,) b. Before 
deriving w.r.t. ky this function, we recall that 



k {n cos P — cos a) + ky (sin a — tan /3 cos a) + O I ky 



Thus, 



Af{ky) 



dk 



y 



sin(Q! — /?) 
cos /3 



b . 



Hence, for the first transmission term, the exit point is 

. , , sin(a — f3) , 

The linear dependence upon b is obvious because of geometry, indeed this expression and those below 
are also derivable by pure geometry. Repeating this procedure for the m-th peak (m = 1 being the 
first), we must include the phase of {R2R1) . This yields 



b] = - 



sin(Q! — f3) 



2 (m — 1) tan /3 cos a 



(15) 



In particular, the resulting values of the peaks for a = 7r/3 arc y^^^^^ = — (2 — m) 5 / \/3- These agree 
very well with the center positions plotted in Fig. 3. The integrated power P of each outgoing laser 
beam relative to the incoming beam power is given by 



2V2 7r 



3/2 



dy 



dky r(0, ky) exp 



^0 ■ ^ \ 



(16) 



For the calculation of the separate laser terms it is sufficient to use the following approxiamtions 
I TjTa I sa 4 n cos a cos (3 j (cos a + n cos /3) and \R2Rx\~\ (cos a — n cos /3) / (cos a + n cos /3) ] . 

These are valid for small k^^y. For example, for the particular case of an incident angle a = tt/3 and 
for a refractive index n = a/S (for which we have a diffusion angle /3 = 7r/6), 

I 1 w 3/4 and | i^^^i | w 1/4 . 

Hence the first three laser transmission beams have the following relative powers 



9 



9 



^.eK'r/3)«(^-,^,^ 

This result is in excellent agreement with the numerical results obtained from Eq.(16) and shown in 
Fig. 3. 



VI. CONCLUSIONS 



Our principal objective in this paper is to develop the formalism, conventions and approximations 
for the interaction of a gaussian laser beam with dielectric structures. We emphasize in particular 
the dominant role played by the momentum perpendicular to the dielectric block, component, 
and in general to any stratified dielectric system. This result draws a strong analogy with one- 
dimensional non-relativistic QM. Indeed we have noted in the text that our system corresponds to a 
QM well potential. The numerical and other calculations displayed in this work invite experimental 
confirmation. It is particularly interesting to note the conditions for coherence and incoherence in 
this study and compare them with those of one-dimensional QM. In QM we deal with time dependent 
wave packets and coherence requires the size of these wave packets to be large w.r.t. the well size. 
In this work the analogy is the separation in the y- variable of laser beams. For this optical study, 
coherence can be achieved not only by increasing the laser beam size w.r.t. the block depth b so that 
overlap occurs between outgoing beams, but also by reducing the incident angle a. This also leads 
to overlap of individual terms since the separation of the y exit points is reduced, for any fixed laser 
size, by decreasing a. We have calculated the various y values for incoherence, both by the SPM and 
numerically and found complete agreement. They remain, of course, to be confirmed experimentally. 
We have also calculated, for various incident angles a, the outgoing incoherent beam intensities and 
profiles. 

A more complex and interesting situation occurs if we consider a dielectric slab cut diagonally, see 
the figure below. 



The dielectric-air-dielectric system now mimics a potential barrier and in addition to diffusion we 
encounter tunneling. Tunneling still holds a number of theoretical conundrums notwithstanding its 
practical importance in technology such as in the tunneling microscope. It seems that technology has 
leap-frogged theory in this case. For example, there is still discussion in the literature about tunneling 
times [23-25] and possible supcrluminal velocities [13,14]. Another theoretical conundrum is that 
for tunneling, our series calculation (two step method) yields a non convergent sum. Transmission 
seems always to be coherent, independent of the barrier/wave packet ratio. In the optical analogy, 
we pass from diffusion to tunneling by simply varying the incident angle. Tunneling occurs when 
<P > /3criticai = arcsin(l/n). 

Let us rapidly recall some of the contents of the previous sections. We started by representing 
the analytic function of a gaussian laser [22] and consequently extended our studies to 3-dimensional 
optics. We also extended this analysis to the laser beam within the dielectric. We then presented 
the optical analogy of a QM well problem and derived both by numerical and approximate analytic 
methods, SPM [20], the exit points from the dielectric block of the laser beam or beams for various 
incidence angles. We have calculated the intensities of the first few incoherent terms (beams) in the 
transition amplitude for a realistic experimental setup. 




Finally, in this section, we have outlined a structure suitable for investigating both diffusion and 
tunneling again dependent upon the incident angle. Another more complex structure would be the 
optical analogy of twin barrier tunneling [16] which, as known, exhibits the extraordinary phenomena 
of resonant tunneling [9,26]. These structures will be the subject of a subsequent report. 
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a = 7r/30 • 


Wgkx '■ 


±2 


±4 


±6 


Wofcy : —6 


(0.345,0.904) 


(0.344,0.904) 


(0.341,0.904) 


(0.337,0.905) 


-4 


(-0.941,-0.284) 


(-0.941,-0.285) 


(-0.940,-0.287) 


(-0.938,-0.291) 


-2 


(0.813,-0.572) 


(0.814,-0.571) 


(0.815,-0.570) 


(0.817,-0.565) 





(-0.056,0.998) 


(-0.056,0.998) 


(-0.059,0.998) 


(-0.063,0.997) 


2 


(-0.746,-0.664) 


(-0.746,-0.664) 


(-0.744,-0.666) 


(-0.742,-0.670) 


4 


(0.978,-0.170) 


(0.978,-0.169) 


(0.979,-0.167) 


(0.980,-0.163) 


6 


(-0.469,0.862) 


(-0.470,0.862) 


(-0.473,0.861) 


(-0.477,0.859) 


q; — tt/6 • 


w„kx '■ 


±2 


±4 


±6 


Wofcj, : —6 


(0.161,0.942) 


(0.160,0.942) 


(0.158,0.943) 


(0.153,0.945) 


-4 


(0.775,0.632) 


(0.774,0.633) 


(0.773,0.634) 


(0.770,0.637) 


-2 


(0.946,-0.024) 


(0.946,-0.023) 


(0.946,-0.020) 


(0.947,-0.016) 





(0.689,-0.518) 


(0.689,-0.518) 


(0.691,-0.516) 


(0.693,-0.511) 


2 


(0.301,-0.764) 


(0.302,-0.764) 


(0.305,-0.763) 


(0.310,-0.761) 


4 


(-0.134,-0.832) 


(-0.134,-0.832) 


(-0.131,-0.833) 


(-0.126,-0.836) 


6 


(-0.624,-0.674) 


(-0.623,-0.675) 


(-0.622,-0.678) 


(-0.629,-0.681) 


a = 7r/3 • 


Wofc:r : 


±2 


±4 


±6 


Wgky : —6 


(0.397,0.480) 


(0.397,0.480) 


(0.395,0.482) 


(0.393,0.484) 


-4 


(0.054,-0.997) 


(0.055,-0.997) 


(0.057,-0.997) 


(0.060,-0.997) 


-2 


(-0.401,0.464) 


(-0.401,0.464) 


(-0.402,0.462) 


(-0.404,0.461) 





(0.670,-0.472) 


(0.671,-0.472) 


(0.671,-0.470) 


(0.672,-0.466) 


2 


(-0.665,-0.345) 


(-0.665,-0.346) 


(-0.664,-0.348) 


(-0.664,-0.350) 


4 


(0.500,0.393) 


(0.500,0.393) 


(0.499,0.395) 


(0.496,0.397) 


6 


(-0.051,-0.983) 


(-0.051,-0.983) 


(-0.049,-0.983) 


(-0.046,-0.984) 



Table 1: Values of the transmission amplitude, (Re[T], Im[r]), as a function of k^ and ky for four 
incident angles, a = (7r/60, 7r/30, 7r/6, 7r/3), for the chosen parameters Wo = 2 mm, A = 633 nm, 
n = and 6 = 5 cm. The dependence of T upon kx is clearly very small. 
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Figure 1: Laser intensity profiles in air for two typical beam waists, Wo = 1 mm and Wq = 2 mm, and 
for a wavelength value of 633 nm. The plot shows the wave number spread and the spatial laser width 
are complementar. 
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Figure 2: Numerical profiles of the modulus of the transmitted (normalized) electric field amplitude, 
for a: = and z = 1 m, as a function of y. Four values of incident angles are shown. For increasing a 
the appearance of secondary peaks is seen. 
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Fig. 3 
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Figure 3: Planar view of the transmitted laser profiles displayed in Fig. 2. The outer curve corresponds 
to an intensity density of 0.01 compared to the intensity of the incident peak, the second curve, where 
it exists, corresponds to 0.05, the third to 0.4 and the fourth to 0.8. For each beam, we also give the 
numerical value of the integrated relative power Prei expressed in Eq. (16). 



